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Abstract 

We present a first-principles study on the dielectric properties of an NaCl crystal and ultrathin 
NaCl films under a finite external electric field. Our results show that the high-frequency dielectric 
constant of the films is not affected by the finite size effect from crystal surfaces and is close to that 
of the crystal, whereas the static one is sensitive to the thickness of the film due to the difference 
in the atomic configurations between the surface and inside of the film. 

PACS numbers: 71.15.-m, 73.61.-r, 77.22.-d 
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I. INTRODUCTION 



The response of insulators and semiconductors to external electric fields is of both fun- 
damental and practical interest. For instance, the response of thin insulator films becomes 
more important as the miniaturization of microelectronics progresses. So far, alkali halide 
has been the subject of a large number of theoretical investigations concerning their dielec- 
tric and piezoelectric characteristics as a typical example of ionic crystals. First-principles 
captions based on the density-functiona. theory QQ have successful* beeu applied to 
studies of the linear response properties of materials at the zero bias limit |3|. However, 
the study of the dynamical transformations of materials in an electric field would greatly 
benefit from the consideration of actual finite fields. Recently, there has been a great deal 
of progress in the first-principles calculation for treating finite electric fields. One method 
is the electronic-structure calculation entirely in real-space, where the boundary conditions 
are not constrained to be periodic. Therefore, one can examine the properties of a material 
by including a finite external field potential in the Kohn-Sham Hamiltonian and imposing 
nonperiodic boundary conditions in the direction parallel to the field Q . Another method is 
the discrete Berry phase scheme , in which the determinations of dielectric 

polarization, Born effective charges, electron-phonon interactions and several other experi- 
mentally measurable quantities of material under a finite external electric field are possible 
using a conventional supercell technique with a three-dimensional periodic boundary. 

In this paper, we implement first-principles calculation to elucidate the dielectric proper- 
ties of an NaCl crystal using the discrete Berry phase scheme and ultrathin NaCl films using 
the supercell with a combination of periodic and nonperiodic boundary conditions. Our 
results reveal that the finite size effect from crystal surfaces on the high-frequency dielectric 
constant of ultrathin NaCl film is negligibly small while the static dielectric constant sensi- 
tively depends on the thickness of the film because of the difference in atomic configurations 
between the surface and inside of the films 

The rest of this paper is organized as follows. In Sec. II, we briefly describe the method 
used in this study. Our results are presented and discussed in Sec. III. We summarize our 
findings in Sec. IV. Finally, mathematical details are described in the Appendix. 
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II. OVERVIEW OF COMPUTATIONAL METHOD 



Our first-principles simulations are based on the real-space finite-difference method 



121 ]. which enables us to determine both the self-consistent electronic ground state and the 



optimised — with a high degree of accuracy by virtue of the t„„e 8 avi„ g double-grid 
technique [13J. The norm-conserving pseudopotentials [14] of Troullier and Martins 



are used as an atomic-ion co re p otential. Exchange-correlation effects are treated by the 
local density approximation [l6| in the den sity functional theory. The nine-point finite- 
difference case, i.e., the case of N=A in Ref [l2j, is adopted for the derivative arising from 
the kinetic-energy operator in the Kohn-Sham equation. A cutoff energy is set at 110 Ry, 
which corresponds to a grid spacing of 0.30 a.u., and a further higher cutoff energy is set at 
987 Ry in the vicinity of the nuclei with the augmentation of double-grid points Q]. The 
structural optimizations are carried out until the remaining forces acting on ions are smaller 
than 33.0 pN. 



III. RESULTS AND DISCUSSION 
A. Dielectric constants of crystal 

The electronic band gap and dielectric constants of an NaCl crystal are calculated. The 
electronic ground state and the optimized geometries under external electric fields are de- 
termined using the discrete Berry phase scheme given by Umari and Pasquarello 9]. The 
calculated and experimental energy differences between the highest occupied and the lowest 
unoccupied levels are listed in the first column of Table [U Since the local density approxi- 
mation underestimates band gaps by approximately 30 %, there is no fundamental reason 
to expect a good agreement between the theoretical and experimental results. The second, 
third, and fourth columns of Table [J show the Born effective charge, and the high-frequency 
and static dielectric constants of the NaCl crystal, respectively. It is well known that when 
an alkali halide is subjected to a static electric field, its polarization becomes larger com- 
pared with the case of the high-frequency electric field because charged ions move along 
the electric field lines. In theoretical calculations, if the ions are kept fixed, this yields an 
electric contribution under an external high-frequency electric field; if both electrons and 
ions are allowed to relax in response to the field, the static dielectric constant is obtained. 
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TABLE I: Calculated and experimental values of the electronic band gap E gap , Born effective 
charge Z*, high-frequency dielectric constant and static dielectric constant eq of NaCl crystal. 





E gap (eV) 


Z* 




£o 


This work 


4.55 


1.11 


2.31 


5.70 


Experiment 


8.6 tt 


1.12 6 


2.25 c 


5.62 c 



"From Ref . 3. 1 • 
^From Ref. E| . 
c From Ref.[ig|. 

Here, to compare the induced charge due to an applied electric field in the crystal with those 
in films in the next subsection, we set the external electric field to 0.25/e^ p V/A, where 
e^ p (=2.25) is the high-frequency dielectric constant obtained in experiment. The electric 
field is applied along the [001] direction of the NaCl crystal. The calculated Born effective 
charge and dielectric constants are in agreement with the experimental values. Figure ^ 
shows the induced charge on the (110) plane due to an applied electric field, in which ions 
are fixed. Since the halogen ions are polarizable units and very little charge remains on the 
Na sites, the primary effect of the applied electric field is the polarization of the CI ions; 
the electrons are pushed from the upper side of the CI ions to the lower side in Fig. ^ and 
across the peak in the 3p charge density of the CI ions. 

B. Dielectric constants of films 

The dielectric constants of ultrathin NaCl films with (001) surfaces are examined. In 
the calculations for the films, the model shown in Fig. El is employed. The supercells are 
taken to be ao x ao x (L + 2)ao and periodic boundary conditions are imposed in the x- 
and y-directions and a nonperiodic boundary condition is imposed in the z-direction, where 
ao(=10.66 a.u.) is the lattice constant of an NaCl crystal obtained in experiments and L 
is the number of lattices in the films. The ionic potential and the ionic Coulomb energy 
under the two-dimensional periodic boundary condition are described in the Appendix. The 
Coulomb charge density cr(= 3.87 x 10~ 4 e/bohr 2 ) on the upper and lower electrodes is 
determined so that the electric field, E, is 0.25 V/A when there is no material between 
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FIG. 1: Induced charge on (110) plane due to applied electric field. Contour interval is 10 
electrons per unit cell. Thick (Thin) curves represent positive (negative) contours and dashed 
curves represent zero charge contours. The atomic sites are indicated for Na by a solid square and 
CI by a solid circle. 

TABLE II: Calculated values of electronic band gap and dielectric constants. 





E gap (eV) 




£o 


L 


= 1 


4.81 


1.95 


2.75 


L 


= 2 


4.69 


2.07 


3.24 


L 


= 3 


4.63 


2.15 


3.62 


L 


= 4 


4.58 


2.21 


3.99 


L 


= 5 


4.55 


2.26 


4.36 



the electrodes. We first examine the optimized atomic geometries of the NaCl films in the 
absence of the external electric field: The Na atom layers on both the surfaces are found 
to sink toward the inside of the films and locate inner than CI atom layers by ~ 0.1 a.u. 
Subsequently, the atomic and electronic structures of the films under the external electric 
field are explored. Figure 01 shows the induced charges on the (110) plane in the cases of 
L — 1, 3 and 5, in which ions are kept frozen. There are no marked differences in the 
polarizations around the CI ions between the thick films and thin films as well as between 
the films and the crystal. 

The computed band gap, high-frequency dielectric constant and static dielectric constant 
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FIG. 2: Computational model for thin films. The film is sandwiched between charged electrodes, 
of the NaCl films are listed in Table |H] The dielectric constants are evaluated according to 

where P and E are the dielectric polarization and the average electric field between the 
electrodes in the ^-direction, respectively [2[J. The band gap of the film becomes smaller 
and approaches that of the crystal, as the thickness of the film increases. The high-frequency 
dielectric constant of the thin films does not significantly vary with the thickness of the films 
and is close to that of the crystal, which account for the results shown in Fig. ^ and Fig. |3] 
where there is no notable difference in the induced charge around the CI ions even though 
the ions locate on the surface of the films. These results represent that the high-frequency 
dielectric constant is dominated only by the polarization of the CI ions. On the other 
hand, the static dielectric constant of the thin films markedly changes according to the film 
thickness. This is interpreted as being due to the difference in atomic configurations between 
the surface and inside of the film; on the surface atomic layer of the upper side, CI ions are 
closer to the upper electrode than the Na ions, while on the other atomic layers, the Na ions 
are closer to the upper electrode than the CI ions (see Fig. 0]). These results indicate that 
the atomic configurations of the films under the finite external electric field play important 
roles in the response of the films to the static electric field. 
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FIG. 3: The induced charge on (110) plane due to the applied electric field for (A) L = 1, (B) 
L = 3, and (C) L = 5. Contour interval is 10 -4 electrons per unit cell. The meanings of the 
symbols are the same as those in Fig. ^ 

IV. CONCLUSION 

We have studied the dielectric properties of the NaCl crystal and ultrathin NaCl films 
under a finite external electric field. When the external electric field is applied, electrons 
are pushed across the peak in the 3p charge density of the CI ions in both the crystal 
and the films. The high-frequency dielectric constant is hardly affected by the finite size 
effects due to the crystal surfaces and is close to that of the crystals. This is consistent 
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FIG. 4: Schematic views of deformations of atomic layers on (010) plane of films for the cases of 
(A) the absence of an external electric field and (B) its presence. Large and small circles are CI 
and Na ions, respectively. 

with the induced charge due to the electric field, in which no significant difference in the 
induced charge around the CI ions is observed between the thick films and thin films as well 
as between the films and the crystal. On the other hand, the static dielectric constant is 
sensitive to the thickness of the film because of the difference in the atomic configurations 
between the surface and inside of the film. The other alkali-halide materials are expected to 
emerge the similar properties to that of NaCl and the work in progress is the investigation 
for them. 
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APPENDIX: COULOMB POTENTIAL AND ENERGY UNDER 2D PERIODIC 
AND ID NONPERIODIC BOUNDARY CONDITIONS 



The ionic potential and ionic Coulomb energy in the cases of two-dimensional periodic and 
one-dimensional nonperiodic boundary conditions are obtained using the Ewald summation. 
Here, the case in which periodic boundary conditions are imposed in the x- and y-directions 
and an isolated boundary condition is imposed in the z-direction is explained. 

By using the identity (a, b, r\ > 0) 



r z e 



2 -a z t 



2+2_ b* 



dt 



'TV 

2b 



e--erfcf^^^ + e-erfc^ + 2 ^ a 



2rj 



2rj 



(2) 



the ionic pseudopotential is given by 



VLir) = ~Z S r^co S [G.(r-R s )]f+(G,r) 



s s 



\ e -\z-R%\W + 0f \ z _ i£| erffl* - R s z \ 77) 



erf(i,|C*|)-^C s ,,erf(^|C ! 



i=l,2 



(3) 



where G = 2tt(^, g-,0), P = (n x L x ,n y L y ,0), S = L x x L y , C = P + r - R s , R s is the 
position of the s-th ion, L x , L y , and L z are the lengths of the unit cell in the x-, y-, and 
^-directions, respectively, and 

/+(G , r)=e -|G|.^. erfc (l G l- 2 ^- fl '') 



+ el G l^'erfc(' G ' +2 y-^). 



(4) 



In addition, C s s {C s -\ +C S)2 = 1) and a s>i are the parameters of the pseudopotential given by 
Hamann et al. [21], the sum of 2~2g^=o ^ s taken over the reciprocal lattice vectors excluding 
the case of \G\ = 0, erf(x) is the error function (or probability integral) defined by erf(x) = 
Jq e~ t2 dt, and erfc(x) = 1 — erf(x). The summations of the first and third terms in Eq. (jSJ) 



are taken over the reciprocal lattice vectors G and the real-space lattice vectors P. When 
r\ is chosen to be 0.2 - 0.7, the number of iterations in Eq. Q is only 7 2 - ll 2 . 
Then, the ionic Coulomb energy is given by 



IE 



2S 



s 



— ^2 Z s Z s t 



1 



\G\ 



cos 



G ■ (R s - R ! 



)} f + (G, R s ') - y~] Z s Z s >5 ss >-^= 
J , V 71 " 



-\R' Z -i*J|V 



V 



erfc(r? |^' s '|) 



(5) 



where $, s ' s = P + R s — R s . A nucleus does not interact with its own Coulomb charge, so 
the P = term can be omitted from the real-space summation when s = s'. The prime in 
the last summation means that |£ s,s | = is omitted. 
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